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ABSTRACT 

Experimental fixed pitch wind turbine performance data is presented 
for both the DOE/NASA Mod-0 and the Danish Gedser wind turbines. 
Furthermore, a method for calculating the output power from large 
fixed pitch wind turbines is presented. Modifications to classical 
blade element-momentum theory are given that improve correlation with 
measured data. Improvement is particularly evident in high winds 
(low tip speed ratios) where aerodynamic stall occurs as the blade 
experiences high angles of attack. 


INTRODUCTION 

Recent tests on the NASA Mod-0 lOOkW wind turbine indicate that 
classical blade element-momentum theory is inadequate when the 
airfoils are at high angles of attack. This problem is particularly 
important in the calculation of fixed pitch and tip control rotor 
performance. Since the maximum power produced by a fixed pitch rotor 
is a critical design parameter, efforts are being made to improve 
theoretical analysis of this operating condition. 

Experimental data from two Mod-0 rotor configurations as well as from 
the Danish Gedser wind turbine were analyzed. An empirical 
correction to the aerodynamic characteristics is presented which 
shows good agreement with the experimental results. 


EXPERIMENTAL DATA 

Tests were run on the NASA Mod-0 100 kW wind turbine to investigate 
the performance characteristics of fixed pitch rotors. Two rotors 
were used in the Mod-0 tests, both of which were 38 meters in 
diameter. Details of the blade planform, twist, etc. are given in 
Tables la and lb and Figures la and lb. Significant differences 
between the two blades are that the aluminum blades h^ve a variable 
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thickness to chord ratio (t/c) and 34 degrees of non-linear twist, 
while the steel spar blades have a constant t/c and no twist. The 
thickness and twist distributions for the aluminum blades are given 
In Figures 2a and 2b. 

The Mod-0 wind turbine was operated with the highly twisted aluminum 
blades at nominal rotor speeds of 20 and 26 rpm. Due to slip of the 
fluid coupling In the drive train, actual rotor speeds were 21.0 and 
27.4 rpm at the maximum power of 50 kW and 105 kW respectively. 
Alternator output power and reference windspeed were recorded on 
magnetic tape and analyzed using the method of bins [ref. 1]. Figure 
3 shows the median measured output power versus reference windspeed. 
The reference windspeed Is measured by an anemometer located at hub 
height about 60 meters upwind of the rotor. This location Is 
believed to give windspeecis which are representative of the average 
freestrearo windspeed at the rotor. The most interesting 
characteristic of the curves Is the relatively constant output power 
at high windspeeds for both rotor speeds. The leveling off of the 
power output occurs for wind speeds at which stall occurs over most 
of the blade. Maximum power Is greater for the higher rotor speed 
since higher wind speed Is required to achieve the same stall angle 
of attack. 

Further tests were conducted using the steel spar blades with no 
twist. The rotor speed was 32 rpm. The inboard 70 percent of the 
blade remains fixed in pitch and thus experiences high angles of 
attack at high windspeeds. As the wind speed increases, the wind 
turbine power output Increases until the generator rating of 100 kW 
Is reached. When wind speed increases further, it Is necessary to 
pitch the outboard blade section so that the generator rating will 
not be exceeded. The tip pitch angle versus nacelle windspeed is 
given in Figure 4. The nacelle wind speed is measured by an 
anemometer located on the wind turbine nacelle. As shown in Figure 
4, the outboard section pitching continues to Increase as wind speed 
Increases. This data indicates a continuing production of torque by 
the inboard fixed pitch portion, of the blade even though that 
section is stalled. 

Finally, performance data from the fixed pitch Gedser wind turbine 
[ref. 2] Is presented in Figure 5. Table 2 and Figure 6 present the 
rotor conf iguratlon for this wind turbine which operated in Denmark 
beginning in 1957. Though not much data was recorded for very high 
winds, operators of the Gedser turbine also reported constant power 
at windspeeds above stall. 

Classical Theory 

Blade element -moment urn theory used in the PROP Code [ref. 3], as well 
as others, divides the blade Into small spanwise elements. These 
elements are each considered to act as airfoil segments in 
two-dimensional flow fields, each at a particular angle of attack. 
From the geometry, the rotational velocity component, the wind 
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component, nnd the "inducinl" oxial and rotational components, the 
local angle of attack is calculated. The lift and drag forces on the 
elements are then determined from two-dimensional ( inf inite aspect 
ratio) wind tunnel data at the local angle of attack fref. 41. 
Comparison of cambered and uncambered data indicates little 
difference in airfoil performance beyond stall. Therefore, data for 
symetric sections may be used [ref. b, hi. The "induced effects" of 
the wind turbine on the flow are determin^'d in an iterative procedure 
until momentum theory is satisfied. The end effects of the finite 
length wind turbine blade are included by using a tip loss model. 
There are a variety of these tip loss models, the simplest being to 
reduce the lift coefficient to zero for approximately 3 percent of 
the radius near the tips. A more complete description of blade 
element-momentum theory, is contained in reference 3. 

The Mod-0 aluminum blade rotor was modelled with the PROP Code. The 
aerodynamic data used is given in Figure 7 for a NACA P3018 
"half-rough" airfoil. The designation "half-rough" denotes 
aerodynamic data which is an average of NACA smooth and rough data 
given in reference 4. This roughness effect was included to account 
for manufacturing imperfections and for the accumulation of dirt as 
the airfoil is exposed to the environment. 

Since the output of the PROP Code is rotor power with no drive train 
losses, the following drive train efficiency model based on 
experimental data was used to calculate the electrical powe>': 


(1) 3 -1.93c + 0.8238 Pj, 

in which P^ is the electrical power, (kW) and P;> is the rotor 
power, (kw). 

Figure 8 shows the predicted power versus windspeed for the Mod-0 
aluminum blades using the PROP Code, The correlation between 
measured data and analysis using infinite aspect ratio airfoil data 
is not very good. Similar results were obtained for the Mod-0 steel 
spar blade [Mgure 9l and the Gedser wind turbine [Figures 10). 

Smooth airfoil data was used for the Gedser wind turbine since 
standard rough data was not readily available for the its airfoil 
(Clark-Y). the most apparent deficiencies of the predicted results 
are (1) the rapid increase in power before stall, {?) the 
less-than-measured maximum power, and (3) the dect'ease in power after 
stall . 



Mode 1 


8y manipulation of the airfoil characterist ics the PROP Code was able 
to match the measured performance of the Mod-0 and Gedser wind 
turbines. A reduction in the unstalled lift curve slope, an increase 
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in the drag coefficient before stall, and a decrease in the drag 
coefficient after stall were required. It was observed that the 
resulting airfoil characteristics were approximately the same as 
those of a finite length wing with the same aspect ratio. 

The formulas for converting infinite length airfoil data to finite 
length data are from the work of Munk, Glauert, and Prandtl. The 
equations are [ref. 7] given by the following: 


(2a) 

Cl 

= c, 
•*0 

(2b) 

Cd 

= Cn + zt(l + a) 
Dq TT/R' 



57.3 C, 

(2c) 

a ‘ 

■ “0 * ■, M 


in which: 

Cl is the lift coefficient 
a is the angle of attack, deg 
Cq is the drag coefficient 
is the aspect ratio 

T a are factors to allow for the change from elliptical span 
loading to an airfoil with rectangular loading 
0 is a subscript denoting infinite aspect ratio data 

The actual aerodynamic load distribution on a wind turbine blade 
varies with windspeed, twist, planform, etc. However, since the 
factors T and a are small, the loading has been assumed elliptical 
( T,a are zero). The above corrections are made to the airfoil data 
below stall. Furthermore, because the end effects of the blade have 
been accounted for in the airfoil characteristics, the tip loss model 
has been eliminated. 

It should be noted that the above corrections are empirical. The use 
of airfoil characteristics for infinite span in classical theory is 
regarded as established. The fact remains, however, that classical 
theory appears incapable of predicting performance at high angles of 
attack which occur for low aspect ratio blades at low rotor tip speed 
ratios. In the extreme, for example, at a tip speed ratio of zero (0 
rpm) classical theory would predict the thrust on the rotor to be 
proportional to the drag coefficient of about 2 for an infinite 
aspect ratio [ref. 5]. We know, however, this is not true. The drag 
coefficient for even a flat plate of aspect ratio 8 is less than 
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1.3 [ref. 8]. It could be that assumption in momentum theory of an 
infinite number of lightly loaded blades needs to be reassessed and 
improved theoretical models developed for this condition. 

Investigation of the airfoil characteristics after stall reveals the 
reason for the predicted negative power at very high winds. Figure 
11 shows an airfoil element operating with its chord line parallel to 
the plane of rotation (as on the steel spar rotor). The resultant 
wind velocity acts at an angle of attack* a , with respect to the 
plane of rotation. The coefficients of lift, Cl, and drag, Cq, 
forces operating on the element can be resolved into a coefficient of 
torque force, Cq , which acts in the plane of rotation. This 
coefficient is given by 

Cq = sin oi - Cq cos a 


Using the infinite aspect ratio data from Figure 7 it can be seen in 
Figure 12 that negative torque can be expected for angles of attack 
between stall and 45 degrees. 


The airfoil characteristics after stall were determined for an 
idealized stall which would result in constant power (torque) at high 
angles of attack. The torque force on a blade element is 
proportional to the coefficient of torque force and the square of the 
resultant velocity or mathematically 


For a constant rotor speed we can divide by the constant which 
yields ^ 


(5) 


q C, 


'R 

4 


but from Figure 11 cosa = / Vj^ 

( 6 ) 

cos^a 

Substituting Eq. (3) gives 


(7) 


C, sin a 

Q ^ - 

cos a 



cos a 


and thus 
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If we let 

( 8 ) 


8 


AjSln 2a + A 


cos^g 
2 sin a 


and 


(9) 


Cq = Bjsln'^a + B2COS o 


and substitute into Eq. (7) we have 


( 10 ) 


Q 'V (2Aj - Bj)s1natan a + (Ag - Bg) 


Since the measured torque after stall Is independent of wind speed it 
must also be independent of angle of attack. Thus, taking the 
derivative with respect to a and setting it equal to zero yields 

(11) Ai»-^ 


Referring to Figure 13, at an angle of attack of 90 degrees equation 
(9) gives 

“ ^DMAX 


For a finite aspect ratio blade 

(13) 1.11 + 0.018 >R 


for ^ < 50 based and experimental data from Reference 8. 


Thus 


(14) 


8 


^DMAX 

2 


Rearranging Eq. (8) and substituting Eq. (14) yields 


(15) 



^OMAX 


sin a cos a) 


sin g 
cos^g 
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Similarly Eqs. (9) and (12) give 


( 16 ) 


Bg = 


' ^DMAX 
cos a 


For continuity with the below stall airfoil data, A 2 and B? 
are solved at the stall angle condition and thus Eqs. (15) and (16) 
become 

sin 

(17) *2 ■ «=Ls' ‘^DMAX “s “s> 


and 


(18) 


In which the subscript s denotes the value of the constant at stall. 

The resulting airfoil characteristics for a NACA 23018 airfoil with 
an aspect ratio of 25 are given in Figure 14 and 15 along with the 
characteristics for an infinite aspect ratio airfoil [Fig. 8]. With 
these corrections to the airfoil characteristics, the predicted 
performance of the Mod-0 aluminum blade is found to correlate well 
with the measured data as shown in Figure 16. Similar results were 
obtained for the Mod-0 steel spar blade and the Gedser wind turbine 
are shown in Figures 17 and 18. 


Bo ® C, 


^DMAX 
cos a_ 


CONCLUSIONS 

A method of correcting the airfoil characteristics for use with blade 
element-momentum theory has been developed. The airfoil data below 
stall is corrected for the finite length of the blade. This approach 
appears to account for the induced effects better than classical 
blade element-momentum theory alone, particularly for highly loaded 
low aspect ratio blades. An idealized model for aerodynamic 
characteristics after stall has been developed which results in 
nearly constant power in high winds. This model shows good agreement 
with experimental data from two rotor configurations on Mod-0 as well 
as the Danish Gedser wind turbine. 
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TABLE la 


ALUMINUM BLADE CHARACTERISTICS 


Rotor d1a,, m(ft) . 
Number of Blades . 
Root cutout, % span 
Fixed Pitch . . . , 

Airfoil 

Taper ....... 

Twist, deg . . . . 

Solidity 

Precone, deg . . , 
Tilt, deg 


38.5 (126.4) 

2 

5 


NACA 230 Series 
.... Linear 
34 (Non-linear) 

.... 0.031 

0 

8.5 


TABLE 1b 

STEEL SPAR BLADE CHARACTERISTICS 


Rotor dia., m(ft) 38.4 (126.0) 

Number of Blades 2 

Root Cutout, % span 23 

Tip Control, % span 30 

Blade Pitch 75% span, deg 0 

Airfoil NACA 23024 

Taper Linear 

Twist, deg 0.033 

Precone, deg 0 

Tilt, deg 8.5 


TABLE 2 

GEDSER BLADE CHARACTERISTICS 


Rotor dia, m(ft) 24 (28.7) 

Number of Blades 3 

Root cutout, % span 25 

Fixed Pitch • — 

Blade Pitch 75% span, deg 6 

Airfoil Clark-Y 

Taper None 

Twist, deg 12 (Linear) 

Solidity 0.090 

Precone, deg 0 

Tilt, deg 10 
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Figure 4 - Measured tip pitch 
angle versus windspeed for 
the Mod-0 steel spar, tip 
control blades 
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Figure 5 - Measured performance 
of the Danish Gedser wind 
turbine 
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Figure 6 - Danish Gedser wind 
turbine blade planform 



Figure 7 - Aerodynamic data for 
"half-rough" NACA 23018 
airfoil 
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Figure 12 - Aerodynamic torque 
force coefficient versus 
angle of attack for "half- 
rough" NACA 23018 airfoil 
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Figure 14 - Comparison of two 
dimensional "half-rough" and 
Improved aerodynamic data 


Figure 13 - General aerodytiamic 
characterist ics for 
idealized stall model 
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Figure 15 - Comparison of two 
dimensional "half-rough" 
and improved aerodynamic 
torque force coefficient 



Figure 17 - Comparison of 

measured and calculated tip 
angle versus windspeed 
using improved aerodynamic 
data for the Mod-0 steel 
spar blades 
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Figure 16 - Comparison of 
measured and calculated 
performance using improved 
aerodynamic data for the 
Mod-0 aluminum blades 



Figure 18 - Comparison of 
measured and calculated 
performance using improved 
aerodynamics for the Danish 
Gedser wind turbines 
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QUESTIONS AND ANSWERS 
L. Vitorna 


From; B. Dohlroth 

Q: Comment on the effect of dynamic stall. 

A: tv's did not ohoewo any oignifioant effeata of dynamic) etall even 

though the 2H0XX aevieo aivfoil were uoed on the MOD-0. Syatem 
dynoinio toada wer>e well behcwcui duning theae tsato. 

From: J. Dugundji 

Q; How did you keep the turbine at constant RPM throughout the wind 
speed range? How do you keep it from overspeeding? 

A; The eleatniaal genenaton ia held at a nean aynohnenoua speed with 
veapeot to the eleotnioal line. The actual notov speed deviatea 
from the nominal RPM due to slip i.. the fluid coupling in the drive 
train. A high speed shaft brake was available to prevent overapeed 
in an emergency upon loss of electrical load. 

From: Anonymous 

Q: If the peak power is more than doubled by going from 20-26 RPM, what 

is the power limit from increasing RPM? 

A: There ia no practical aerodynamic limit to maximum power with in- 

creasing RPM. At higher rotor speeds the rotor stalls at a higher 
windspeed and thus a much higher power. 

From; M. Rolland 

Q: How would you characterize the starting ability of the fixed pitch 

rotor? Is there a sacrifice with fixed pitch? 

A: This ia not addressed in the presentation. We have studied starting 

cnaracteristics of fixed pitch rotors and believe it is^ possible to 
design a fixed pitch, machine which will start in low winds and not 
be penalized significantly in performance. The zero twist blades 
on MOD-O start with about 4 degrees of pitch in a 6 m/s wind. Ibe 
S4 degree twisted blades start with 0 degrees of pitch in a similar 
wind. 

From: K. Foreman 

Q: How do you explain the revised aero-characteristics of the blade 

after stall, also Reynolds number effects? 

A: We believe this is due to three-dimensional flow effects which are 

not accounted for using the classical blade element-momentum theory. 
The airfoil data is modeled at the Reynolds number at the 75 per- 
cent span of the blade. 
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Vltorna tt^ontlnuod) 


From: (!. R. Kotloy 

Q; Do you have uny oxporloncc of tho effort on the stalloU power curve 
of the smooth surface- laminar flow aerofoils, instead of the half- 
rough 23XXX series? 


A: }\n<j'oma>hU' J\'i> if NAVA tiivj\uU /a 

f'i'aa//<t itf'i’ iiiiU'H i-u Aohn a’/iiaiZi'W'e "iUaU-bh{Ui.‘(\l hint cihil- 

! tij oj' a '/Vi'/i'f'i'i/ Roti'v. " 
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